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Q' 

A complete description of Q-conditional symmetries for two classes of reaction- 

. diffusion-convection equations with power diffusivities is derived. It is shown that all the 

\ known results for reaction-diffusion equations with power diffusivities follow as particular 

' cases from those obtained here but not vise versa. 

> ■ 

\ 1. Introduction. 

o . 

^ ■ Nonlinear reaction-diffusion-convection (RDC) equations of the form 

^ : 

O; Ut = [A{U)U^]^ + B{U)U, + C{U), (1) 

Oh' where U = U{t, x) is the unknown function, A(U), B(U), C{U) are given smooth functions and 

r-| ■ the subscripts t and x denote differentiation with respect to these variables, arise in a wide range 
, of mathematical models describing various processes in physics and biology [U El [3] . Starting 
^ \ from the remarkable Ovsiannikov work |i4] a great number of papers devoted to investigation 
\ of these equations by means of group-theoretical methods. At the present time one can claim 
that all Lie symmetries of ([1]) are completely described and the relevant Lie solutions are 
constructed for equations of the form ([1]), which arise in applications (see [5]-[T0] and the 
papers cited therein). 

In 1969 Bluman and Cole [TT] introduced an essential generalization of Lie symmetry using 
the simplest representative of ([T]), the linear heat equation. These generalized symmetries are 
often called nonclassical symmetries nevertheless this notion was not used in [11]. The notion 
of nonclassical symmetry was further developed in [12] and [13]. A new generalization of Lie 
symmetry, conditional symmetry, was suggested by Fushchych and his collaborators [H], [13 
Section 5.7]. Note that notion nonclassical symmetry can be derived as a particular case from 
conditional symmetry but not vise versa (see, e.g., an example in [IS]). In the middle of 90-s 
of the last century the notion of generalized conditional symmetry was introduced [17], which 
again can be considered as a special case of conditional symmetry. Taking this into account, 
to avoid any misunderstanding we continuously use the terminology Q-conditional symmetry 
instead of nonclassical symmetry. In fact, there are several types non-Lie symmetries at the 
present time and each of them can be called nonclassical one. 

While there is no existing general theory for integrating nonlinear RDC of the form ([1]), 
construction of particular exact solutions for these equations is a non-trivial and important 
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problem. Finding exact solutions that have a physical, chemical or biological interpretation is 
of fundamental importance. It is well-known that the notion of Q- conditional symmetry plays 
an important role in investigation of nonlinear RDC equations since, having such symmetries 
in the explicit form, one may construct new exact solutions, which are not obtainable by the 
classical Lie machinery. Several papers were devoted to this topic during the last 15 years 
P, [13 [ISl [ISl EDI ED [22] • The time is therefore ripe for a complete description of non-Lie 
symmetries for the general RDC equation ([T]). Since it seems to be extremely difficult task at 
the present time here we present the solving for some important particular cases of ([1]), namely: 

Ut= [U^U^l + XU'^U^ + CiU), (2) 

Ut = [U^U^l + Af/^+^t/, + C{U), (3) 

where A and m are arbitrary constants while C{U) is an arbitrary functions. 

It should be noted that Q- conditional symmetry of ([2]) with A = 0, m = 0, i.e. reaction- 
diffusion equation, was investigated in [15], [18], the most general results were obtained in [201 
[23] . Operators of Q-conditional symmetry for equation (j2]) with A = 0, m have been 
constructed in [21] while the complete description of the RDC with m = is presented in 
the recently published paper [22]. Finally, we remind the reader that the determining equations 
for constructing the Q-conditional symmetry operators of the general RDC equation ([1]) were 
obtained in [9j, however, that paper contains only examples of particular solutions of those 
equations. 

We stress that only reaction-diffusion equations with the convective terms (A ^ 0) are 
considered below. The motivation of this restriction has two aspects. The first one is to 
find Q- conditional symmetries for nonlinear equations involving three transport mechanisms 
(diffusion, reaction and convection) in contrary to standard reaction-diffusion (RD) equations. 
The second one is to deal with the equations, which arise in applications. In fact, sometimes 
the convection arises as a natural extension of a conservation low and then one obtains RDC 
equations instead of RD equations. The effect of nonlinear convection in RD equations can 
have "a dramatic effect on solutions" [21 Section 11.4]). 

The paper is organized as follows. In the second section, we present two theorems giving 
a complete description of Q-conditional symmetries of the nonlinear RDC equations ([2])-(l3]). 
In the third section, the proof of the theorems are presented. The main results of the paper 
are summarized and discussed in the last section. In the second part of this paper, we shall 
apply the Q-conditional symmetries obtained for constructing new exact solutions of the RDC 
equations arising in various applications. 

2. Main results. 

We want to find all possible Q- conditional symmetries of the form 

Q= dt + at,x,U)d, + r]{t,x,U)du, (4) 

where ^ and r] are unknown functions, for the RDC equations ([2])-(l3]). We don't consider the 
problem of constructing Q-conditional symmetries of the form 

Q = da: + ri{t,x, U)du, 
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because one is equivalent (up to the known non-local transformation) to solving the given 
equations ©-([3]) [23]. 

Now we present main results of the paper in the form of two theorems. Note that we 
search for purely conditional symmetry operators, which cannot be reduced to Lie symmetry 
operators. 

Theorem 1 Equation ^ is Q- conditional invariant under the operator if and only if it and 
the relevant operator (up to equivalent representations generated by multiplying on the arbitrary 
smooth function M{t, x, U) ) have the following forms: 

{i) Ut = [U^Ux]. + XU'^Ux + (Aif/'^+i + A2)(f/-"^ - As), m ^ -1, A2 ^ (5) 

Q= dt + {XiU + X2U-'^)du] (6) 

(n) Ut= [U-'Uxl + XU-'Ux + {XihiU + X2){U - Xs), Xiy^O, (7) 
Q= dt + (XilnU + X2)Udu; (8) 

(m) Ut= [U-"Wxl + XU~'^Ux + XiU + X2U-^ + A3, (9) 

Q= dt +f{t,x)dx +2{g{t,x)U + h{t,x)u'^)du, (10) 
where the function triplet (f,g,h) is the general solution of the system 

Vf. + ft + fg = o, 

fxx - Xfx - 2gx - fh = 0, 

i9-^)ig + 2fx)+gt = 0, (11) 
2gh - Xih + 2fxh - X2fx + ht - Xqx - Qxx = 0, 
h^-^h- X^fx + ^g- Xhx - hxx = 0. 

Hereafter Ai, A2 and A3 are arbitrary constants. 

It should be noted that cases (i) and (ii) with A = immediately give the RD equations 
and the relevant symmetries obtained in [2T] . 

Nevertheless system ( ITTll contains five equations on three unknown functions it is com- 
patible. In fact, the system with f = g = 0, Ai = is reduced to the ordinary differential 
equation 

hxx + Xhx + y /l - /l' = (12) 

and hence 

Q= dt +2h{x)U^du (13) 



is the Q-conditional symmetry operator for an arbitrary non-zero solution of fll2p . Unfortu- 
nately, ODE (1121) cannot be integrated for the arbitrary coefficients A and A2, however, some 
particular solutions can be easily established. For example, setting h = a particular case of 
(i) with m = — |, Ai = is obtained. 

Setting A = A2 = in f|T2l) . we arrive at the known ODE hxx = h"^ with the general solution 
h = yV(0, Ci, X-I-C2), where Ci and C2 are arbitrary constants, W is the Weierstrass function with 
the periods and Ci. Its simplest solution takes the form h = 6x~^ and leads to the known 
Q-conditional symmetry operator Q = dt + 12x~'^U2du of the nonlinear diffusion equation 
Ut = [U^2Ux]x [21]. However the result derived in [21] can be generalized as follows. One can 
easily check that an arbitrary particular solution of (fTT!) with A = generates the Q-conditional 
symmetry operator fllOl) of the RD equation 

Ut= [U-'Wx]x + XiU + X2U-^ +A3. (14) 

Obviously, the operator presented in Table 3 of [2T] is obtainable from (ITU]) and (ITT!) by setting 
A = Ai = and f = g = but not wise versa. 

Theorem 2 Equation ^ is Q-conditional invariant under the operator ^ if and only if it and 
the relevant operator (up to equivalent representations generated by multiplying on the arbitrary 



smooth function M{t, x, U) ) have the following forms: 

(t) Ut = [U^'UxU + XU'^^'Ux + XiU + m + -1, (15) 

Q= dt-XU"'+'dx + {XiU + X2U~"')du; (16) 

(n) Ut= [U-'^Uxl + XU'^Ux + iX,u'^ + X2U-^+Xs)Q^ + U^y (17) 

Q= dt+ (^-A[/^ + ^) + (^1^' + ^2^' + (18) 

(m) Ut = Uxx + XUUx, (19) 

Q = a, + Qf/ + dx+(^a + bU~ ^U' - ^U'^ du, (20) 

where triplet of the functions (a, b, q) is the general solution of the system 
o-t = CLxx ~ 2ag^, 

bt = bxx - 2bqx + Aa^, (21) 
qt = Qxx - '2qqx - 26^; 

(zv) Ut = Uxx + XUUx + Xo + X2U\ A2^0, (22) 



4 



Q = dt 



-xu 



A 



a. + (Ao + \2UM 



•u, 



(23) 



(v) Ut = U,. + XUU, + Xo + XiU + XsU\ Aay^O; (24) 

Q^ = dt+p^Ud^ + -^^{Xo + XiU + X3U^)du, 2 = 1,2 (25) 

where pi are the roots of the quadratic equation 2p^ + Xp + 9A3 — A^ = 0, 

Q = dt + bd, + (7fo.. - KU)du, (26) 
where the function b{t, x) is the general solution of the overdetermined system 
(A7 - 3)b,, + 2bb^ + bt = 0, 

bxxx - bbxx + Xibx = ^, (27) 
X-ibxxx + bl + 3Xibx + ^b = Co, 

where 7 = g^,AoAiA3 7^ and Aq, Cq G M. 

Remark 1. The case of the Burgers equation (fT9|) was completely investigated in paper 
[25] (see also a locally equivalent equation to the Burgers equation found in [H] and particular 
solutions of presented in [HISS]), while the cases {iv) — (v) were obtained in [22] • 

Remark 2. The RDC equations 
Ut = Uxx + XUUx + Ao + Alt/ + A2/72 

and 

Ut = Uxx + XUUx + Ao + Alt/ + Ast/' + X3U' 

possess Q-conditional symmetry, however they are reduced to ([22]) and (^^, respectively, by 
the simple local substitutions [22]- 

Inserting A = into flM|) - fl25l) one immediately arrives at the well-known RD equation and 
Q-conditional symmetry operator with cubic nonlinearities constructed earlier [T8]-|20]. 

We have also proved that the first two equations of fl271) with A = and fef = is equivalent 
to the equation 66^. — 26^ + 9Ai = derived in ^20] (see P. 264) to construct the Q-conditional 
symmetry operator fl2Bl) in the explicit form for the RD equation with cubic nonlinearity. 

Consider some equations which arise as particular cases of those from Theorem 1 and 2 and 
are known in application. Equation ([5]) with m = 1 contains as a subcase the equation 

Ut= [UUx]x + XUUx + U{l-U), (28) 

which is a natural generalization of the equation 

Ut= Uxx + XUUx + U{1-U), (29) 



extensively studied by Murray Section 11.4]. On the other hand, fl28p is nothing else but 
the porous-Fisher equation with the Burgers convective term UUx- Note that the Murray 
equation fl29|) arises in Theorem 2 (see case (iv) and Remark 2). Its Q-conditional symmetry 
was established earlier in [9] and several exact solutions were recently found in [22]. Equation 
([9]) with A2 = -Ai = -2, A3 = 0, t 2t takes the form 

Ut= (t/^U + A(/7^), + f/5(l-t/5). (30) 

This equation may be called the Murray equation with the fast diffusion. Another analog of 
the Murray equation with the fast diffusion is 

Ut= (t/^U + A[/^t/, + t/5(l-f/5), (31) 

which is a particular case of equation f|T5l) . 

Consider the generalized Fitzhugh-Nagumo (FN) equation 

Ut = U,. + XUU, + XsUiU - 6){1 -U), 0<5<1. (32) 

We remind the reader that (152]) with A = is the famous FN equation [27| describing nerve 
impulse propagation. It can be also considered as a simplification of Hodgkin-Huxley model 
(see, e.g., [3]) describing the ionic current flows for axonal membranes. Equation fl32]) with A = 
and 5 = 1 is the Kolmogorov-Petrovskii-Piskunov equation, which firstly was investigated in 
[28] (see English translation in [29|) to describe the population dynamics under some restrictions 
on characteristic individuals. Equation (1521) can be reduced to the form 

Wt = Wyy + XWWy + Aq + A 1 1^ " X^W^", (33) 

where 

Ai = A3(^(5+ 1)^-5), Ao = A3^(5+l)(^(5 + 1)^-5), (34) 

by the local substitution 

W{t,y) = U-^{5+l), y = x + ^{5+l)t. (35) 

Now one note that equation (1551) is nothing else but equation (|25p with the new notation. On 
the other hand, (ITT)) contains as a particular case the equation 

Ut = + A/7^f/. + U"^m - 6){l -U), 6= ^, (36) 

which may be treated as a generalized FN equation with the fast diffusion. In the case A = ±3, 
this equation may be called the generalized Kolmogorov-Petrovskii-Piskunov equation with the 
fast diffusion. 



3. Proofs of the theorems. 

Proof of Theorem 1. The proof of Theorem 1 and 2 is based on the known algorithm 
for finding Q-conditional symmetry operators (see, e.g., [9], [I5]). Firstly, we apply the local 
substitution 

f ?7™+\ m ^ -1, 

V = \ ' ^ ' (37) 
[ In [/, m = —1. 

In the cases m ^ —1 and m = — 1, substitution (1571) reduces equation ([2]) to the forms 
V^^ = V''Vt-XV^ + F{V), (38) 

(here n = 0, F{V) = -{m + l)C(K^), A 7^ 0) and 

V,, = exp{V)Vt-XV, + F{V), F{V) = C{expV), (39) 

respectively. 

The determining equations for the general RDC equation 

V,, = Fo{V)Vt + Fi{V)V, + F2{V), 

being Fi(y),i = 1,2,3 arbitrary functions, have been obtained in ^ (see P.535). In the case 
Fq{V) = y", Fi{V) = -A and F2iV) = F{V) those equations take the form 

Cvv = 0, 

= 2ey(-A-e^") + 2e.y, 
{2^vV - - 6)V^" - CvnV^-' - A^,. + 3^^^ - 2rj,v + U = 0, 
r]Fv + (2e. - r]v)F + nrj^V"-' + 2^,r]V^ + r]tV'' - Xr], - r],^ = 0. 

Solving the first equation of fHUl) . we arrive at the function ^ = a{t,x)V + f{t,x) being 
a{t, x) and f{t, x) arbitrary smooth functions at the moment. 

It turns out that system (140|) doesn't possess any Q-conditional symmetry if a{t,x) ^ 0. So 
we must assume 

^ = f{t,x). (41) 
Solving the second equation of (HUi) under condition (HTl) . we arrive at 

T] = g{t,x)V + h{t,x). (42) 
Taking into account fHTj) and (j^2l) the third equation of fHOl) reduces to the form 

(2//.. + ft + nfg)^^ + nfhV'^-' - + A/,. + 2^,. = 0. (43) 

This equation can be splitted with respect to the powers of V. One needs to consider two cases 
depending on n: 
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(a) if n 7^ 1 then 

2//. + ft + nfg = 0, 

fh = 0, (44) 

fxx — ^fx ~ '^9x = 0. 

(b) if n = 1 then 

Vfx + ft + fg = o, 

fxx - Xfx - 2g^ -fh = 0. 

Let us consider case (a). Substituting (HTj) and (142|) into the fourth equation of (140|) one 
arrives at 

(gV + h)Fv + (2/, -g)F = -nV^'\gV + hf + + Xh,+ 

H9xx + Xgx)V - {gt + 2f,g)V^+^ - {K + Vxh)V^. ^ ' 

To solve (jlSD and (jUj) one needs to consider two subcases, which follow from the second equation 
of dm), i.e. either / = or /i = 0. 
The case / = leads to the system 

/ = 0, 

9x = 0, (46) 
{gV + h)Fv -gF = -nV^-\gV + hf + + \K - gtV^^^ - h^V^. 

Setting g = const, h = const, we arrive at the system 

f = o,g = xi,h = x;,F = (xiv + a;)(A3 - r"), 

therefore 

V,. = V^'Vt - AK + {XIV + A;)(A3 - Vn, (47) 

Q = dt + {XIV + X;)dv. (48) 

Applying substitution (!37|) to the equation (H71) and operator psj) we obtain case (z) of the 
theorem (note one should use new notations Aj = i = 1,2). 

Now we assume that g ^ constant, so that the third equation of (l46l) can be reduced to the 
form 



V + -] Fv-F = -ngV''-' {V+-] + i _ _ (49) 

9/ \ 9/ 9 9 9 

It turns out that the last equation can be satisfied only under condition ^ = const (see the 
proof below). Setting ^ = constant into fH9l) and making the relevant calculations, we obtain 
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only Lie symmetry operators and a particular case of operator fl48p and equation (H7|) . For 
example, if ^ = then the system 

h = 0, 

2ffx + ft + nfg = 0, 

fxx ^fx '^dx Oi 

gVFv + (2/,. - g)F = -nV^+^g^ + + Xgx)V - {g^ + 2f,g)V^+^ 



(50) 



is obtained, which can be easily solved and its general solution has the form 

/ = ^'7'^'"') . , = „ = 0, F = A.l/»« + A.V, 

C2 exp[Aint) + 1 C2 exp(Aini:j + 1 

where G M, k = 1,2. Hence we arrive at the RDC equation 

Vxx = V^Vt - XVx + XiV^+' + X2V 
and the operator 

Ciexp(Aint) AiV" 

Q = ^t + 7V—^V^^^ T^^V^^y- 

C2 exp(Ainrj + 1 C2exp(Aim) + l 

However one can establish by multiplying (15T|) on the function M{t,x,U) = 1 + C2exp(Ai?2t) 
that the last operator is nothing else but the Lie symmetry operator (see case 8 of Table 1 in 

Let us prove that ^ = const. By differentiating equation fl49p with respect to the variables 
X and t one obtains two equations. Assuming {^)t{^)x = 0, one easily arrives at the condition 

- = const. 

9 



' h\ ( h 

X we arrive at the equation 



Consider the case {-)t{-)x 7^ 0. By differentiating equation (H9|) with respect to the variables 



1/^ , .t\Trn ^TT ^n—l ^xxx ~(" Xhx 



-{2ngh + ht)xV'' - 2nhV''-' + 



X X 

Since the function V"', V"^~^ and 1 on the right-hand-site are functionally independent (we 
consider the case n 7^ 1) their coefficients must by constants. It means that 2nh = const so 
that hx = 0. Taking now two last equations of fj44l) . one easily establish that gx = 0, i.e. we 
arrive at the contradiction: (^)j. = 0. 

Consider the case (b). Substituting (HTj) and fH2]) into the fourth equation of (HU]) . we arrive 
at (1451) with n = 1. Dealing with this equation in the same way as above (see the case (a)) we 
obtain equation 

Vxx = vvt - XVx + xiv^ + x;v + x;, (52) 

and the operator 

Q = dt + fit, x)dx + {git, x)V + h{t, x))dv, (53) 
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were the triplet (/, g, h) are the general solution of fllip . Applying formula fl37p with m ^ —1 
we obtain the case [iii) of the theorem (note one should use new notations Aj = — 2A*, i = 
1,2,3). 

Finally, we analyze equation fl39|) . which is locally equivalent to the RDC ([2]) with m = —1. 
Using again the determining equations, which have been obtained in |9] to find operators of the 
Q-conditional symmetries (P), we arrive at the following system 

Cvv = 0, 

Vvv = 2^y(-A - ^exp V) + 2^^y, 

(6 + 2ee. - 2CvV + ^V) exp V + XC, - 3^vF + 2r],,v - U = 0, ^^^^ 
rjFv + {2^x - Vv)F + {rf + 2^^?] + r]t) exp V - Xiq^ - rj^^ = 0, 

where ^, rj and F are yet-to-be determined functions. Solving the first and second equations of 
this system we establish that the functions ^ and r] must be given by formulas (HTIl and (H2|) . 
respectively, otherwise Q-conditional symmetry doesn't exist. Substituting fHT]) and fH2l) into 
the third equation of we obtain the equation 

ifh + ft + 2ff,) exp V + {fg)Vexp V + Xf, + 2g, - f,, = 0. 

Since the functions /, g and h don't depend on V, one can split this equation with respect to 
exp V and V exp V and obtain the system 

fh + ft + 2fU = 0, 

fa = 0, (55) 
A/.. + 2g, - = 0. 

Substituting fj4T|) and fj42l) into the fourth equation of fl54|) we arrive at the equation 

(gV + h)Fv + (2/,. - g)F = -{gV + hf exp V + /i,.,. + A/i,+ 
+ (^7xx + A(?,.)^ - + 2/.^)!^ exp V-{ht + 2f,h) exp V. ^^''^ 

Now we apply to (IH^ the same approach, which has been used for solving equation 0451) . Thus 
taking into account system 0551) . we obtain finally the expressions 

F= (AiV + A2)(A3-expy), / = 0, ^ = Ai, h = X2, 

which lead to the equation 

= exp{V)Vt - AK + [XiV + A2)(A3 - exp V) (57) 

and the operator 

Q = dt + {X,V + X2)dv. (58) 

Applying substitution (1371) with m = —1 to (1571) and (158|) one obtains the case (n) of the 
Theorem 1. 
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The proof is now completed. 

Proof of Theorem 2. First of aU we note that all Q-conditional symmetries of equation 
([3]) with m = were found in the recent paper [22] so that the restriction m 7^ is assumed 
below. We again use substitution fl37|) . which reduces equation ([3]) to the form 



= V'Vt - Ay"+V, + FiV), (59) 
(here n = ^ 0, -1, F{y) = -(m + l)C(V^)) if m ^ -1 and to the form 

V;, = exp{V)Vt - Xexp{V)V, + F{V), F{V) = C(exp V) (60) 
if m = —1. 

Consider equation (!59|) . Using the general form of the determining equations obtained in 
[9j one easily arrives at the following system 

Cvv = 0, 

Vvv = '2Cv{-XV^+'-CV^)+2C.v, 

r]Fv + (2e. - r]v)F + nrfV^~^ + 2^.77!/- + r/^V^" - \V^+^r^^ - r],, = 0, (61) 
Aex.^"+'+((-2^y + A(n + l))r] + 2^^, + ^t)^^ + ^r]nV''-^ - 3^vF+ 

to find the function ^, t] and F. In the case of equation ( I6OII that system takes the form 
Cvv = 0, 

r]vv = -2^1/ (A + exp V + 2^.^, 

(6 + 2ax. + (A + e - 2ey)^ + A^,.) exp V - S^yF + 2r/,y - = 0, ^^^^ 
77Fy + (2^^ - r]v)F + {if + 2^^?] + rjt - Xr]^) exp F - 77^^ = 0. 

Taking into account the first equation in fl6Tl) and fl62|) . one establish that there are only 
three possibilities for the functions ^ and r] : 

{a)c = \iv + x;, v = viv), xi,x;em, 

(b) e = fit, x), ri = g{t, x)V + hit, x), (63) 

(c) ^ = a(t, x)V" + fit, x), T] = r]it, X, V), ait, x) 7^ 0. 

In the case (c) function 'qit, x, V) takes the forms 

- (;S|S3) ^"""^ - (Trrfife) ^"''^ + + 9it, x)V + hit, x),n^ -2, -3, 
V = \ -2a(a + A)\/ In y + 2a/ In y + a^V^ + (2a(a + A) + git, x))V + hit, x), n = -2, (64) 
2a(a + A) In y - afV'^ + a^V^ + git, x)V + hit, x), n = -3; 
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for the system (IMi) and 



-2aVexp V - 2a{\ + f -2a) exp V + a^V^ + g{t, x)V + 



for the system 

Consider case (a) and system flM]) . Since the function ^ is the hnear function, the general 
solution of the second equation of (16T!) is the function 



(n + 2)(n + 3)' (n + l)(n + 2)' ' ' ' ^^^^ 

where A3, A4 G M, n ^ —2, —3. Substituting (165!) into the fourth equation of (1611) . one obtains 



3At 



Ai(n - 2) + A(n + 1) H/" + nA^^ 



1 *T /n— 1 



(66) 



if A^ 7^ (the case A* = leads only to the Lie symmetry operators). Substituting fl^Bl) into 
the third equation of system flMl) one arrives at the expression 



(1 + ^(AI(n - 2) + \{n + 1))^"'^ + (n - 1)A;F 



*T An— 2 



0. 



The first possibility is 



1 



1 + ^(^K^ - 2) + Kn + 1) ) V"-' + {n- 1)A;\/"-2 







(67) 



(6^ 



while the second one is = 0. 

Splitting (168|) with respect to the different powers of V one obtains the system 

(AJ + A)(n + l) = 0, 
(n-l)A^ = 0. 



(69) 



The first equation immediately gives A]; = — A because n 7^ — 1. The second equation of fl69l) 
can be transformed into an equality if A2 = or n = 1. Hence we obtain the equation 



= V^Vt - Ay"+V, - {KV + \l)V\ 
and the relevant Q-conditional symmetry operator 

Q = dt- XVd, + (A*V^ + \l)dv 
if A2 = 0, and the equation 

= vvt - xVv. - ^i±|^^ ( \\i\v'^ + \iv + a: 



3A 



and the relevant operator 



Q = 9i + (-Ay + A;)a,. + ( -\\\v' + A*y + a: ) dy, 



(70) 
(71) 

(72) 
(73) 
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if n = 1. Finally, applying substitution (^71) with ra ^ —1 to (170]) - (175]) and introducing new 
notations A3 = Ai(m + 1), A4 = A2(m + 1) in the case of equation (170]) and Ag = A3 = 
^7 A4 = ^ in the case equation (172]) . we arrive exactly at the items {%) and {ii) of the Theorem 
2. 

It turns out that the special values n = — 2 and n = — 3 (see ( 164]) ) lead only to particular 
cases of equation ( 170]) and operator ( 17T]) . Moreover the case = again leads to the equation 
(170]) and operator ^ with A^ = A^ = 0. 

To complete the investigation of case (a) one needs to consider system (162]) . Integrating the 
second equation one easily obtains 

ri = 2A*(-A*'^ + 2A* - A - A^) exp V + X^V + A*, (74) 

where AJ, A2, A3, A4 G M. Inserting ^ = XIV + X2 and ([71]) with A^ 7^ into the third equation 
of ([62]), we find 

F = ^^^^(>^l^ + >^ + >^*2-^K). Kt^O. (75) 
Substituting ( 175]) into the fourth equation of ( l62]) . we arrive at the condition 

+ 3I* + A + a; - At)^ exp V = 0, (76) 

which can be satisfied only in the case i] = 0. However it leads to the requirement A^ = (see 
( 174]) ) what contradicts to ( 175]) . 

Finally, system ( 162]) with A^ = takes the form 



(77) 



^ = xiv + xi, 

iX*,V + Xl)iX*, + X)=0, 
{X^V + Xl)Fv - A^F = -{X*^V + XlY exp V. 

Solving (j77]) with respect to the function F, one obtains the equation 

V;, = exp(T/)T4 - Aexp(\/)\4 + (X^V + Xl){X; - exp V) (78) 
and the relevant Q-conditional symmetry operator 

Q = dt-Xd^ + {X;V + Xl)dv. (79) 

Applying substitution (137]) with m = — 1 to (178]) and (179]) and introducing new notations 
A3 = Ai, A4 = A2, A5 = A3, we obtain the RDC equation 

Ut = [U~'U4, + XU^ + (Ai \nU + X2){U - A3) (80) 

and the Q-conditional operator 

Q= dt-Xd, + iXi\nU + X2)Udu. (81) 

1.-^ 



However equation (ISUj) and operator (ISTj) are reduced to the reaction-diffusion equation 

Ut = [U-'Uy]y + (XilnU + X2)iU - X3) (82) 
and the operator 

Q= dt + {\i\nU + \2)Udu (83) 

by the local substitution y = x + Xt. Note that fl82l) and fl83l) were previously derived in [2T] . 

Consider case (b). It turns out that this case leads only to Lie symmetry operators. Indeed 
the first and the second equations of system fl^ are automatically satisfied. Substituting flB51) 
into the fourth equation of (I6T1) . one obtains 

A ((n + l)g + /,) (A(n + l)/i + /, + 2//, + 72/(7) + ri/W^-^ + 2(7,. - = 0.(84) 

Now one needs to consider two subcases: n ^ 1 and n = 1 (we remind also the reader that 
n ^ —1,0). Assuming n 1 and splitting (18^ with respect to the different powers of V, we 
obtain the system 

{n+l)g + f,, = 0, 

X{n + l)h + ft + 2ff, + nfg = 0, 

fh = 0, 

"^Qx fxx 0- 

The third equation of fl85p leads to / = or /i = 0. Setting / = one immediately arrives at 
g = h = since n + 1 and A 7^ 0. Hence we obtain a particular case (a). Setting h = and 
using the first and the fourth equations of fl85|) we obtain 

9 = 9{t), f = -{n + l)g{t)x + ^{t), 

where the functions g{t) and Lp{t) must be determined. Substituting g and / into the second 
equation of (ISBj) and splitting with respect to x, the ODE system 

gt-in + 2)g' = 0, 
(ft - {n + 2)Lpg = 

with the general solution 

9 = 7 ^ , f = 7 ^ , Cfc G M, A; = 1, 2 (86) 

^ (n + 2)t + ci' ^ (n + 2)t + ci' ' ' ^ ' 

is obtained. To find the function F{V) we solve the third equation of fl6T|) . taking into account 
the expressions found for /, g and h, and obtain F = XiV'^"''^^. Hence we arrive at the equation 

= V'Vt - Ay"+V, + Aiy2"+^ (87) 



u 



and the operator 



g = g,+ l" + ^^" + ^^ 4-7 ^ dy. (88) 

^ {n + 2)t + ci {n + 2)t + ci ^ ' 

However the operator (l88l) is not purely Q-conditional one since it is equivalent to the operator 

D=(^{n + 2)t + ci) dt+ ((n + l)x + C2) - Vdy. (89) 

Applying substitution fl37|) with m 7^ — 1 to fl87|) and (l89l) one arrives at the equation and a 
linear combination of Lie symmetry operators listed in the case 7 of Table 1 [S] . 
Assuming n = 1, we obtain only a particular case of (1571) and fISS]) . 

Solving system (!62l) we again find only Lie symmetry operators of equation fl60l) with 
F = Ai + A2 exp and F = Ai exp(A2l^). 

Dealing in a quite similar way with the case (c) we have established that this case doesn't 
produce any new Q-conditional operators. 

The proof of Theorem 2 is now completed. 



4. Conclusions. 



In this paper, Theorems 1 and 2 giving a complete description of Q-conditional symmetries 
of the nonlinear RDC equations ([2]) - ([3]) are proved. It should be stressed that all Q-conditional 
symmetry operators listed in Theorems 1-2 contains the same nonlinearities with respect to the 
dependent variable U as the relevant RDC equations. Analogous results were earher obtained 
for single reaction-diffusion equations [15], [H], [19], [20] . 

However, we note that there is the essential difference between RDC equations ([2]) - ([3]) and 
the relevant RD equation 

Ut = [U"'U,l + CiU). 

For example, the Murray type equation fl22|) admits the Q-conditional symmetry fl23l) . while 
the RD equation with this term, i.e. the Fisher type equation 

Ut = U,. + Ao + XiU + A2f/^ A2 ^ 

does not possess one. Similarly, the RDC equation (fT7|) possessing the Q-conditional symmetry 
ffTSjl doesn't has an analog among reaction-diffusion equations with the diffusivity f/~2. 

The RDC equations listed in Theorems 1 and 2 contain several well-known equations arising 
in applications and their direct generalizations. In the particular case, the Murray equation 
fl2U]) . its porous analog (1251) and its analogs fl5Ul) - (l^ with the fast diffusion; the Fitzhugh- 
Nagumo equation [27] with the convective term 

Ut = a,^, + XUU, + XsUiU -5)il-U), 0<S<1, (90) 

and its analog ( |36l) with the fast diffusion; the Kolmogorov-Petrovskii-Piskunov equation [28] 
with the convective term 

Ut = u,. + xua, + x^u{i - uf, (91) 



and the Newell- Whitehead equation [30] with the convective term 

Ut = U,, + XUU^ + X^U' - X,U. (92) 

A further generalization of the RDC equations ([2]) and reads as 

Ut = [U^U,l + XU^U, + C{U), A^O, (93) 

where m and n are arbitrary constants. It turns out that all possible Q-conditional symmetries 
of fl93l) coincide with those presented in Theorems 1 and 2. In other words, if a nonlinear 
equation of the form (l93l) admits a Q-conditional symmetry operator ( not Lie symmetry !) 
then either n = m, or n = m + 1. The work is in progress on the complete description of 
Q-conditional symmetry of the RDC equation 

Ut = [exp{mU)U^]^ + Xexp{nU)U^ + C{U), X^O. 

It is well-known that new Q-conditional symmetries don't guarantee the construction of 
exact solutions, which cannot be obtained by the Lie machinery (see non-trivial examples in 
[311 [32]). Ill the second part of this work we will demonstrate that the Q- conditional symmetries 
obtained above can be successfully applied for constructing new non-Lie solutions. Moreover we 
will show that some of those solutions have remarkable properties. In the particular case, it will 
be shown that those solutions may satisfy typical boundary conditions arising in mathematical 
biology. 
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